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Abstract  Matched  field  processing  (MFP)  provides  a  means  of  attaining  the  full 
gains  available  from  the  shallow-water  acoustic  channel  in  passive  sonar  signal 
processing.  By  modeling  the  full  field  structure  of  acoustic  signals  propagating  in 
the  ocean,  MFP  offers  the  potential  for  both  detection  gain  (through  its  better  signal 
model)  and  localization  gain  (through  its  additional  discrimination  capability  in  range 
and  depth)  over  traditional  planewave  processing.  However,  high  spatial  ambiguities 
and  mismatch  present  formidable  challenges  in  practice,  limiting  the  performance 
gains  that  are  realistically  achievable  with  MFP. 

Prediction  of  MFP  localization  performance  is  a  challenging  problem.  MFP  replica 
(steering)  vectors  can  be  highly  ambiguous  in  range  and  depth,  resulting  in 
significant  non-local  estimation  errors  at  low  signal-to-noise  ratios  (SNRs) — errors 
not  modeled  by  traditional  localization  measures  such  as  the  Cramer-Rao  bound. 
Recent  work  has  demonstrated  the  accuracy  of  an  interval-error-based  method, 
referred  to  herein  as  the  "method  of  interval  errors"  (MIE),  in  predicting  mean- 
squared  error  localization  performance  well  into  the  threshold  region  where 
non-local  errors  may  dominate. 

This  work  uses  the  MIE  to  predict  the  mean-squared  error  accuracy  of  MFP  range 
and  depth  estimates  for  two  well-known  approaches:  (i)  conventional  beamforming 
(equivalent  to  maximum  likelihood  estimation  for  white  noise)  and  (ii)  Capon-MVDR 
adaptive  beamforming.  Simulation  results  will  characterize  localization  performance 
as  a  function  of  SNR,  for  apertures  and  environments  of  interest.  Particular  attention 
will  be  given  to  the  "threshold  SNR"  (below  which  localization  performance  degrades 
rapidly  due  to  global  estimation  errors)  and  to  the  minimum  SNR  required  to  achieve 
acceptable  range/depth  localization.  Initial  work  will  also  be  presented  assessing 
the  MIE's  potential  to  characterize  localization  performance  in  the  presence  of 
mismatch. 
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ABSTRACT 

The  full-field  modeling  of  matched  field  processing  (MFP) 
provides  potentially  large  gains  over  traditional  beamform- 
ers  in  passive  sonar  signal  processing.  MFP  localization 
gains  are  manifested  by  superior  estimation  of  source  range 
and  novel  estimation  of  source  depth.  However,  MFP  local¬ 
ization  is  limited  in  practice  by  high  spatial  ambiguities  in 
the  MFP  output.  These  ambiguities  can  result  in  global  lo¬ 
calization  errors  at  low  signal-to-noise  ratios,  errors  which 
must  be  accounted  for  to  predict  MFP  localization  perfor¬ 
mance  accurately.  This  work  uses  the  so-called  method  of 
interval  errors  (MIE)  to  predict  MFP  mean-squared  error 
localization  performance  well  into  the  “threshold  region” 
where  global  errors  dominate.  New  results  enable  highly 
accurate  threshold  region  MIE  predictions  for  both  conven¬ 
tional  beamforming  and  Capon-MVDR  adaptive  beamform¬ 
ing  in  the  simple  signal-in-white-noise  case.  A  variation 
of  the  standard  MIE  technique  is  shown  to  produce  equally 
accurate  performance  predictions  when  both  mismatch  and 
colored  noise  are  introduced. 

1.  INTRODUCTION 

Matched  field  processing  (MFP)  provides  a  means  of  attain¬ 
ing  the  full  gains  available  from  the  shallow- water  acous¬ 
tic  channel  in  passive  sonar  signal  processing.  By  model¬ 
ing  the  full  field  structure  of  acoustic  signals  propagating 
in  the  ocean,  MFP  offers  the  potential  for  both  detection 
gain  (through  its  better  signal  model)  and  localization  gain 
(through  its  superior  estimation  of  source  range  and  novel 
estimation  of  source  depth)  over  traditional  planewave  and 
range- focused  beamforming. 

MFP  localization  performance  is  limited  in  practice  by 
high  redundancy  in  the  MFP  steering  vectors,  which  results 
in  spatial  ambiguities  in  MFP  output.  These  ambiguities, 
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FI 9628-00 -C-0002.  Opnions,  interpretations,  conclusions,  and  recommen¬ 
dations  are  those  of  the  authors  and  are  not  necessarily  endorsed  by  the  U.S. 
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which  are  often  signficant  in  the  range  and  depth  dimen¬ 
sions,  can  cause  large  localization  errors  at  low  signal-to- 
noise  ratios  (SNR’s).  Localization  is  defined  here  as  esti¬ 
mation  of  a  source  location  parameter  (e.g.,  source  range, 
depth,  or  bearing)  by  maximizing  beamformer  output  over 
the  parameter  search  space.  Localization  error  is  defined  as 
the  physical  distance  between  the  estimated  parameter  value 
and  the  true  parameter  value,  quantified  here  in  a  mean- 

squared  sense:  E  |  ^0  —  Otruej  |  ,  where  9  is  the  estimate 

of  the  true  parameter  9 true  >  and  where  E  is  expectation. 

For  example,  consider  the  MFP  range- depth  beampat- 
tern  shown  in  Fig.  1,  for  source  location  6  km  range  and 
25  m  depth.  At  high  SNR’s,  range/depth  estimates  for  the 
source  will  likely  fall  within  the  “local  error”  region  sur¬ 
rounding  mainlobe  of  the  beampattem.  As  SNR  decreases, 
however,  sidelobe  ambiguities  become  more  likley  to  cause 
global  errors ,  where  the  source  is  localized  outside  the  beam- 
pattern  mainlobe.  Several  potential  global  error  regions  are 
marked  in  the  Fig.  1. 
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Figure  1:  MFP  range-depth  beampattem  for  a  horizontal 
line  array  (length  600  m)  at  endfire,  A/2  spacing  at  fre¬ 
quency  50  Hz.  Environment  is  off  the  southern  California 
coast.  True  source  location  is  6  km  range  and  25  m  depth. 
Grey  scale  units  are  dB  (relative  to  peak). 


Mean-squared  error  performance  curves  are  typically  char¬ 
acterized  by  three  distinct  (SNR)  ’’regions”  of  performance, 
as  depicted  in  Fig.  2:  the  asymptotic  region  in  which  lo¬ 
cal  errors  dominate  and  MSE  is  proportional  to  log(SNR); 
the  threshold  region ,  in  which  global  errors  dominate  and 
MSE  rises  rapidly  as  SNR  decreases;  and  the  no  informa¬ 
tion  region ,  in  which  the  signal  is  below  the  noise  and  lo¬ 
cation  estimates  are  essentially  uniformly  distributed.  Two 
important  SNR’s  of  interest  are  the  threshold  SNR ,  below 
which  global  errors  become  significant,  and  the  required 
SNR ,  the  SNR  that  is  needed  to  achieve  minimum  accept¬ 
able  MSE  performance.  These  SNR’s  are  important  for  sys¬ 
tem  design,  as  will  be  described  below.  For  matched  field 
processing,  acceptable  performance  boundaries  often  occur 
within  the  threshold  region,  meaning  that  accurate  thresh¬ 
old  region  performance  prediction  is  vital  to  understanding 
what  SNR’s  are  truly  required  for  acceptable  MFP  localiza¬ 
tion.  Ignoring  global  errors  for  a  local  error  bound  such  as 


Figure  2:  Typical  mean-squared  error  (MSE)  performance 
vs.  phone-level  SNR.  Asymptotic  region  is  all  SNR’s  above 
-5  dB.  Threshold  region  is  between  -18  and  -5  dB.  No  in¬ 
formation  region  is  below  -18  dB. 

the  Cramer-Rao  bound  leads  to  falsely  optimistic  estimates 
of  the  SNR  required  to  achieve  acceptable  performance.  As 
can  be  seen  in  Fig.  2,  if  the  linear  part  of  the  MSE  curve  is 
extended  all  the  way  across  the  figure,  the  intersection  with 
the  acceptable  performance  line  would  occur  around  -40  dB 
SNR,  much  lower  than  the  more  accurate  estimate  of  -10  dB 
required  SNR  that  takes  global  errors  into  account. 

Clearly,  then,  quantifying  MFP  localization  performance 
is  a  challenging  problem.  What  is  needed  is  accurate  thresh¬ 
old  region  MSE  performance  prediction  that  takes  into  ac¬ 
count  all  of  the  following:  ambiguities  in  the  MFP  range - 
depth  beampattern  (that  lead  to  global  errors  at  low  SNR’s), 
signal  model  mismatch,  colored  noise,  and  finite-sample  ef¬ 
fects  (for  adaptive  beamformers).  Previous  methods  that 
account  for  global  errors  have  involved  calculating  bounds 
which  give  only  approximate  performance  predictions  in  the 
threshold  region  [2,  71  or  have  focused  on  (non-adaptive) 


maximum  likelihood  estimation  in  white  noise  with  known 
data  covariances  [1,  8].  The  approach  taken  here  is  to  ap¬ 
ply  new  results  [4]  to  an  established  technique  [5],  referred 
to  herein  as  the  method  of  interval  errors  (MIE),  to  en¬ 
able  accurate  threshold  region  MFP  performance  prediction 
for  two  popular  beamformers  in  passive  sonar,  the  adaptive 
Capon-MVDR  (CAP)  beamfonner  and  the  non-adaptive  con¬ 
ventional  beamformer  (CBF).  Section  2  describes  MIE  and 
the  new  results  that  enable  its  use  for  CAP  and  CBF.  Sec¬ 
tion  3  presents  some  simulation  results  demonstrating  the 
accuracy  of  this  approach  for  threshold  region  performance 
prediction  and  also  describes  a  new  procedure  to  extend 
the  MIE  framework  to  include  both  signal  model  mismatch 
and  colored  noise.  Section  4  concludes  by  summarizing  the 
findings  of  this  work. 

2.  APPROACH 

2.1.  Framework 

This  work  assumes  a  stochastic  signal  model,  as  is  com¬ 
monly  adopted  in  passive  sonar.  Data  snapshots  x  are  as¬ 
sumed  to  be  zero-mean,  complex-Gaussian:  x  ~  Cr\  (0,  R)  . 
The  covariance  R  =  cr^vsv^  +  cj^Rn  ,  where  v5  is  the  true 
signal  vector,  is  the  signal  power,  cr^  is  the  noise  power, 
and  Rn  is  the  noise  covariance,  normalized  so  that  its  trace 
equals  |vs|2  and  cr^jc is  the  phone-level  SNR.  The  sam¬ 
ple  covariance  matrix  R  is  then  formed  in  the  usual  way: 

L 

R  =  (l/L)  E  XfcXj[*  ,  where  L  is  the  number  of  snapshots. 
k=i 

This  analysis  focuses  on  estimators  that  are  derived  from 
two  popular  beamformers  in  passive  sonar,  the  adaptive  Capon- 
MVDR  (CAP)  beamfonner  and  the  non-adaptive  conven¬ 
tional  beamfonner  (CBF).  Each  estimate  is  obtained  by  max¬ 
imizing  the  beamfonner  power  output  over  values  of  the  pa¬ 
rameter  6  in  some  search  region  0: 


#CBF 

=  argmaxCBF(0) 

=  argmaxwgBF(e)RwCBF(0) 

(i) 

OcAP 

=  argmaxCAP(0) 

ote 

=  argmaxwgAPWRwCAP(0) 

=  argmax  (vf  R-1v4  , 

0<E©  V  / 

(2) 

where  the  CBF  weight  vector  Wcbf(0)  =  /  | |2  is  sim¬ 

ply  a  scaled  version  of  the  steering  vector  v  0  at  6  and  the 
CAP  weight  vector  is  a  function  of  the  sample  covariance: 

WCAP(0)  =  R.-1v<?/  (vf  R_1V^  .  (3) 

The  maximum  likelihood  estimate  for  the  stochastic  sig¬ 
nal  model  described  above  may  be  derived  in  the  case  where 


the  underlying  covariance  R  is  known  but  the  SNR  is  un¬ 
known.  The  MLE  is  given  by 

Omle  =  arg  max  [MF  (0)  -  log  (MF  (0))]  ,  (4) 

where  the  matched  filter  output  wmf(0)-^wmf(0)  *s  a  func¬ 
tion  of  the  matched  filter  weight 

wMf(0)  =  .  (5) 

Note  that  (5)  depends  on  knowledge  of  the  underlying  noise 
covariance  Rn  ,  which  cannot  be  estimated  in  passive  sonar 
because  signal-free  training  data  is  never  available.  Thus, 
the  MLE  serves  mainly  as  a  benchmark  for  MFP  perfor¬ 
mance  in  this  paper.  Note,  however,  that  the  MLE  is  equiv¬ 
alent  to  CBF  in  the  case  of  white  noise  (Rn  =  I)  and  known 
covariance. 


2.2.  Method  of  Interval  Errors 


The  method  of  interval  errors  decomposes  the  expression 
for  mean-squared  error  into  two  conditional  probabilities, 
that  an  interval  error  (IE)  has  occurred  or  that  no  interval  er¬ 
ror  (NIE)  has  occurred,  “interval  error”  being  synonymous 
with  “global  error”  : 


£  j(<?-0tr«e)2} 

=  Pr  (NIE)  E I  (e  -  etrue )2  NIE  j 
+Pr(IE)E  j  (O-Otme)2  IeJ 

Nm 

1  -  E  P  [9 = M 

m= 0 

H"  ^  ^  P  \  @  =  @mj  X  (9m  9true) 


x  MSE, 


asym 


(6) 


m— 0 


where  MSEas?/m  represents  the  asymptotic  (high-SNR)  MSE 
and  the  set  1  represent  parameter  values  correspond- 

ing  to  Nm  side  lobe  peaks  in  the  beampattern.  As  seen  in  Eq. 
(6),  MIE  approximates  the  probability  of  a  global  error  by 
recognizing  that  in  the  threshold  region,  global  errors  will 
almost  always  occur  at  sidelobe  peaks  in  the  beampattern. 
This  approximation,  which  turns  out  to  be  highly  accurate 
under  a  wide  variety  of  conditions,  turns  what  would  be  a 
continuous  integral  into  a  finite  sum.  A  further  savings  is 
attained  by  employing  a  union  bound  approximation: 


p(o  =  Kpe=p  (BF  (0m)  >  BF  ( Otrue ))  ,  (7) 


where  BF  ( 9 )  is  the  beamformer  power  at  position  9  (e.g., 
CAP  (9)  or  CBF  (0)),  and  where  the  right-hand  side  of  (7) 


is  the  dominant  pairwise  error  probability  pe  in  the  union 
bound.  Substituting  (7)  into  (6)  gives  the  form  of  MIE  used 
here. 

Using  MIE  thus  requires  two  algorithm-specific  quan¬ 
tities:  the  asymptotic  MSE,  MSEas?/m ,  and  the  pairwise 
error  probability  pe.  The  expressions  for  MSEasym  are  pre¬ 
viously  established  results.  For  MLE,  MSEasym  is  given 
by  the  familiar  Cramer-Rao  bound  (CRB),  which  may  be 
computed  using  the  Slepian-Bangs  formula  and  depends  on 
R-1  and  the  vector  partial  derivative  with  respect  to  pa¬ 
rameter  9 ,  d\s/d9  .  For  CAP  and  CBF,  the  expressions  for 
MSEasym  can  be  derived  using  Taylor  series  expansions 
of  the  beamformer  output  power  about  the  source  location 
9 true  =  9 s  and  are  given  in  works  by  Vaidyanathan  and 
Buckley  (VB)  [6]  and  Hawkes  and  Nehorai  (HN)  [3],  re¬ 
spectively.  Both  expressions  depend  on  the  data  R,  the  first 
derivative  dvs/d9 ,  and  the  second  derivative  d2vs/d9 2  , 
and  the  VB  expression  also  depends  on  the  third  deriva¬ 
tive  d3vs/d9 3  .  Because  MFP  steering  vectors  are  super¬ 
positions  of  normal  mode  functions  that  are  environment- 
dependent,  the  vector  derivatives  cannot  be  represented  an¬ 
alytically  and  must  be  computed  numerically  using  finite 
differences.  For  the  results  presented  in  this  paper,  it  was 
found  that  a  step  size  of  0.01  m  was  sufficient  to  compute 
the  numerical  vector  derivatives  accurately  for  both  range 
and  depth.  (Note  that  the  derivatives  need  only  be  computed 
at  the  source  location.) 

The  important  results  that  allow  application  of  MIE  to 
both  CAP  and  CBF  are  the  pairwise  error  probabilities  pe 
recently  derived  by  Richmond  [4]  for  CAP  and  CBF.  Again, 
p€  represents  the  probability  that  the  beamformer  power  at 
some  location  9m  will  exceed  that  power  at  the  true  source 
location  9 true  =  9S.  The  procedure  for  computing  pe  for 
CAP  is  as  follows  (the  derivation  is  found  in  [4]): 

1.  Define  the  TV  x  2  matrix  V  =  [vm|vs]  ,  where  vm  is 
the  steering  vector  corrsponding  to  the  sidelobe  peak 
9m  ,  then  compute  the  QR -decomposition  of  R-1/2  V  , 

D2x2 

0(AT-2)x2  J 

and  let  di  and  d2  denote  the  two  2x1  columns  of  D. 

2.  Compute  the  eigenvalues  of  the  2x2  matrix  d2d^  — 
did^  and  denote  the  eigenvalues  as  Ai  and  A2. 

3.  The  error  probability  pe  for  CAP  is  then  given  by 


pe  =  0.5  •  [1  +  sgn  (Ai)] 

— sgn  (Ai.)  •  T (— A2/Ai.,  L  —  N  +  2)  , 


r-i/2v  _  Q 


Mean-Squared  Error  for  Depth 


is  the  cumulative  distribution  function  for  a  special 
case  of  the  complex  central  F  statistic. 

To  compute  pe  for  CBF,  follow  the  same  procedure  as  above 
with  the  following  modifications:  replace  R-1/2  by  R1/2 
in  Step  1,  use  the  2x2  matrix  did^  —  d2d^  (instead  of 
d2d^  —  didf*)  to  compute  the  eigenvalues  Ai  and  A2  in 
Step  2,  and  replace  L  —  N  +  2  by  L  in  Step  3. 

3.  SIMULATION  RESULTS 

The  simulation  results  that  follow  are  for  a  horizontal  line 
array  with  AT  =  41  elements  and  15  m  spacing  (total  aper¬ 
ture  600  m).  The  array  was  assumed  straight  and  its  depth 
was  assumed  to  be  100  m.  MFP  replicas  were  generated 
using  the  KRAKEN  normal  mode  model  for  an  environ¬ 
ment  near  the  Santa  Barbara  channel,  with  depth  568  m 
and  downward -refracting  sound  speed  profile.  Beamformer 
outputs  (from  which  MFP  range  and  depth  estimates  were 
derived)  were  computed  over  a  search  interval  of  2-10  km 
range  and  1-101  m  depth,  at  endhre.  Beamformer  outputs 
were  generated  using  a  sample  covariance  R  with  L  =  60 
snapshots  (L  ss  1.5 IV). 


Figure  3:  MFP  mean-squared  error  for  range  estimation 
versus  phone -level  SNR  for  signal  in  white  noise  (no  mis¬ 
match)  scenario. 

In  the  first  simulation,  a  50  Hz  endhre  source  was  placed 
at  6  km  range  and  25  m  depth  and  the  noise  was  assumed  to 
be  white  (Rn  =  I)  and  Gaussian  with  power  cr2  =  70  dB. 
In  the  white  noise  case,  parameter  estimates  derived  from 
CBF  and  from  MLE  are  identical.  Figs.  3  and  4  show  the 
prediction  results  for  estimation  of  MFP  range  and  depth, 
respectively.1  The  results  compare  the  MIE  predictions  of 

1  In  this  paper,  range  and  depth  were  estimated  separately,  with  the  other 
parameter  assumed  known.  MLE  has  a  straightforward  extension  to  simul- 


Figure  4:  MFP  mean-squared  error  for  depth  estimation  ver¬ 
sus  SNR  for  signal  hi  white  noise  (no  mismatch)  scenario. 

MSE  for  CBF  and  CAP  (as  a  function  of  phone -level  SNR) 
to  MSE  curves  generated  from  4000  Monte  Carlo  runs  at 
each  SNR  point.  MIE  predictions  for  both  range  and  depth 
and  both  CBF  and  CAP  are  accurate  for  all  SNR’s,  but  most 
importantly  the  predictions  are  accurate  within  1-2  dB  in 
the  threshold  region.  Note  that  in  the  absence  of  interferers, 
CBF/MLE  outperforms  CAP  because  CAP  suffers  losses 
due  to  finite-sample  adaptive  training.  Boundaries  on  ac¬ 
ceptable  values  of  \/MSE  were  taken  nominally  to  be  500 
m  in  range  (i.e.,  significantly  better  than  the  mainlobe  range 
resolution  of  range-focused  beamforming  at  the  same  fre¬ 
quency,  «  2.5  km  at  2  km  range)  and  10  m  in  depth  (i.e., 
enough  to  enable  discrimination  of  surface  and  submerged 
sources).  Note  again  that  these  boundaries  fall  well  within 
the  threshold  region  for  both  range  and  depth.  Thus,  pre¬ 
dicting  the  threshold  region  MSE  well  is  vital  to  predicting 
the  required  phone-level  SNR  (denoted  SNRr  here)  for  ac¬ 
ceptable  MFP  localization.  In  Fig.  3,  SNRr  is  «  —14.5  dB 
for  CBF  and  &  — 11  dB  for  CAP,  while  in  Fig.  4,  SNRr  is 
«  -17.8  dB  for  CBF  and  «  -14.8  dB  for  CAP. 

Because  MIE  predictions  are  based  on  the  beampattem 
of  a  source  at  a  specific  location,  it  is  important  to  under¬ 
stand  how  performance  varies  as  source  location  varies.  To 
do  this,  one  may  re-compute  the  curves  of  Figs.  3  and  4 
for  each  source  range  and  depth  in  the  search  space,  derive 
SNRr  for  both  range  and  depth,  and  use  the  higher  of  the  two 
SNRr  values  as  the  overall  SNRr  at  a  given  source  location. 
For  CBF  in  the  current  scenario,  this  procedure  reveals  a  3- 
5  dB  variation  in  SNRr  over  the  range -depth  search  space. 
From  the  values  of  SNRr,  one  can  back  out  the  source  lev¬ 
els  required  for  acceptable  MFP  localization  by  adding  back 


taneous  estimation  of  range  and  depth,  but  simulating  the  latter  is  much 
more  time-consuming. 


the  noise  level  of  70  dB  and  the  predicted  transmission  loss 
(TL)  at  each  source  location.2  This  computation  is  shown  in 
Fig.  5  for  CBF.  Note  that  the  source  levels  shown  in  Fig.  5 
are  reflective  of  the  ideal  nature  of  this  scenario,  white  noise 
and  no  mismatch. 
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Figure  5:  Required  source  levels  for  acceptable  MFP  local¬ 
ization  using  CBF,  as  function  of  source  range  and  depth. 
Gray  scale  units  are  dB  re:  1  ^Pa/Hz. 

When  signal  model  mismatch  is  introduced  into  the  prob¬ 
lem,  the  MIE  procedure  needs  to  be  modified  to  account  for 
the  possibility  that  9  may  never  converge  even  asymptoti¬ 
cally  (with  increasing  SNR)  to  9 true  [4].  Most  often,  9  will 
converge  to  some  mismatched  value  9  mm  and  there  will 
be  an  asymptotic  bias  in  MSE  equal  to  1 9 true  ~  Omm  f  ■ 
To  account  for  this,  one  needs  to  compute  MSEasym  us¬ 
ing  the  true  data  R  and  the  mismatched  vector  v  mm  (cor¬ 
responding  to  location  9mm)-  Further,  the  necessary  error 
probabilities  pe  describing  the  transition  to  the  threshold  re¬ 
gion  are  now  derived  from  the  competition  of  the  local  side- 
lobe  peaks  of  the  “mismatched  beampattem”  (computed  by 
taking  normalized  inner  products  of  the  mismatched  steer¬ 
ing  vector  with  the  true  source  vector)  with  the  peak  mis¬ 
matched  vector  v mm  ■ 

Thus,  the  modified  MIE  expression  for  MSE  is  given 
by  (6)  and  (7),  except  with  9  mm  replacing  9 true  in  the  pair- 
wise  error  probabilities  and  with  MSEasym  computed  using 
vmm  instead  of  the  true  source  vector  vs.  The  interval  er¬ 
ror  is  still  computed  relative  to  9 true  in  order  to  reflect  the 
actual  accounting  of  MSE.  And  the  error  probabilities  are 
still  computed  using  the  procedure  detailed  in  the  previous 
section,  with  vmm  replacing  vs  in  Step  1. 

hi  the  second  simulation,  a  mismatch  in  compressional 
sound  speed  of  10  m/s  was  introduced  into  the  propagation 
model  (true  value  was  1572,  mismatched  value  was  1562); 
for  this  case,  the  noise  was  still  white.  Fig.  6  shows  the 
prediction  results  for  estimation  of  MFP  range  in  the  mis- 

2The  MFP  propagation  model  gives  TL  for  each  phone  in  the  array, 
and  one  can  use  the  root-mean-squared  of  all  the  phone  TL’s  to  derive  an 
overall  ’’array  TL”  at  a  given  source  location. 


matched  case  (depth  estimation  performance  was  very  sim¬ 
ilar).  Note  that  the  mismatch  causes  a  slight  degradation 
in  performance  in  the  threshold  region  and  an  asymptotic 
bias  of  «  4  m  (5  dB).  Importantly,  however,  the  threshold 
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Figure  6:  MFP  mean-squared  error  for  range  estimatiion 
versus  SNR  for  mismatched  signal  in  white  noise  scenario 
(mismatch  of  10  m/s  in  compressional  sound  speed). 

region  predictions  are  still  accurate  within  1-2  dB  for  both 
CBF  and  CAP.  Because  this  scenario  again  involved  white 
noise,  CBF  still  outperformed  CAP  because  of  the  adaptive 
training  losses  incurred  by  CAP. 

Finally,  one  can  apply  the  modified  MIE  procedure  for 
the  mismatch  case  to  the  colored  noise  (discrete  interferer) 
scenario,  for  cases  where  Rn  is  fixed  (but  not  white)  and 
where  all  discrete  interferers  are  located  outside  the  search 
region.3  In  the  specific  problem  examined  here,  Rn  = 
cr^dd^  +  cr2 1 ,  where  d  represents  the  steering  vector  of  a 
distant  interferer  and  where  the  interferer  power  a  \  is  held 
constant  while  the  signal  power  cr2  (and  thus  the  signal-to- 
ambient-noise  ratio  fSANR])  varies.  Constant  power  for 
distant  inteferers  is  not  unusual  for  passive  sonar  data,  so 
this  is  a  case  of  interest. 

In  this  colored  noise  case,  MSE  performance  for  very 
high  SANR’s  will  minor  the  asymptotic  region  of  the  white 
noise  case,  because  the  signal  is  dominant  and  parameter 
estimates  will  be  near  the  signal  location.  As  the  signal 
power  cr2  (and  thus  the  SANR)  decreases,  the  signal-to- 
interference  ratio  cr2  /cr|  will  decrease  and  the  interferer  will 
become  more  influential,  with  increasing  likelihood  that  pa¬ 
rameter  estimates  will  be  near  the  largest  sidelobe  of  the 
interferer  within  the  search  space  (recall  that  the  interferer 
itself,  d,  is  located  outside  the  search  space).  Denote  the 

3If  interferers  are  located  within  the  search  region,  then  the  problem 
becomes  a  multi-source  parameter  estimation  problem.  MIE  with  slight 
modification  still  produces  accurate  performance  predictions  in  this  case, 
as  detailed  in  [4]. 


steering  vector  associated  with  this  largest  interfere!*  side- 
lobe  as  Vrf.  Then  the  analogous  threshold  region  MSE  per¬ 
formance  can  be  derived  by  following  the  same  procedure 
as  in  the  mismatch  case  outlined  above,  with  replacing 
vMm  ■  The  added  complication  here  is  that  the  SANR  de¬ 
pendence  of  the  “colored  beampattem”  (computed  by  tak¬ 
ing  nonnalized  output  of  the  beamformer  weights  with  the 
data,  both  computed  using  ideal  covariances)  needs  to  be 
accounted  for. 

For  the  third  scenario,  an  interferer  was  introduced  at 
20  km  range,  5  m  depth,  and  10  0  bearing,  with  power  g\  — 
66.8  dB  (corresponding  to  150  dB  source  level).  Fig.  7 
shows  the  prediction  results  for  estimation  of  MFP  range 
in  the  colored  case  (depth  estimation  performance  was  very 
similar).  Note  that  with  the  introduction  of  the  interferer, 
CBF  is  no  longer  equivalent  to  MFE,  and  CBF  performance 
is  degraded  significantly,  by  about  10  dB  SNR  relative  to  the 
white  noise  case.  CBF  still  outperforms  CAP  in  the  asymp¬ 
totic  region  (where  behavior  is  like  the  white  noise  case), 
but  the  adaptive  CAP  now  outperforms  CBF  in  the  thresh¬ 
old  region,  where  adaptive  interference  rejection  becomes 
more  beneficial.  Note  that  threshold  region  performance  for 
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Figure  7:  MFP  mean- squared  error  for  range  estimatiion 
versus  SNR  for  signal  in  colored  noise  scenario. 

CAP  is  still  well-predicted  by  MIE  within  1-2  dB,  while 
CBF  tines  hold  region  prediction  is  off  by  slightly  more  be¬ 
cause  significant  interference  side  lobes  cause  greater  depar¬ 
tures  from  MIE  assumptions.4 

4.  SUMMARY  AND  DISCUSSION 

This  work  demonstrates  a  powerful  prediction  tool  for  MFP 
localization  performance.  Mean-squared  error  (MSE)  pre- 

4Threshold  region  performance  for  MLE  is  also  approximate  because 
precises  error  probabilities  were  not  computed. 


diction  was  shown  to  be  accurate  within  1-2  dB  in  the  thresh¬ 
old  region  for  both  MFP  range  and  depth  and  for  both  the 
CBF  and  Capon-MVDR  (CAP)  beamformers.  Accurate  thresh¬ 
old  region  predictions  in  turn  lead  to  better  estimates  of 
the  SNR’s  (and,  equivalently,  the  source  levels)  required  for 
acceptable  MFP  localization.  These  accurate  threshold  re¬ 
gion  predictions  for  both  CBF  and  CAP  were  made  possible 
by  application  of  new  pairwise  error  probability  results  to 
the  established  method  of  interval  errors  (MIE)  technique. 
Slight  modification  of  the  basic  MIE  procedure  also  enables 
accurate  threshold  region  prediction  in  the  cases  of  signal 
model  mismatch  and  colored  noise,  as  was  demonstrated  in 
separate  simulation  examples. 

It  is  important  to  note  that  because  MIE  is  situation- 
specific  (for  example,  one  needs  to  hypothesize  a  mismatch 
mechanism  for  the  signal  model  mismatch  case),  defini¬ 
tive  quantification  of  MFP  localization  performance  will  re¬ 
quire  systematic  analysis  repeating  the  above  simulations 
over  several  reasonable  mismatch  and  interference  scenar¬ 
ios.  Future  work  may  also  extend  the  techniques  presented 
here  to  the  case  where  sources  are  moving  during  the  obser¬ 
vation  interval. 
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The  MFP  Localization  Problem  4g£ATO 


Shallow-Water  Multipath  Propagation: 


models  acoustic  multipath  propagation 
to  enable  3-D  source  localization 


•  Goal:  quantify  what  SNR’s  are  required 
for  acceptable  MFP  localization 

•  Approach:  accurate  MFP  mean-squared 
error  prediction  that  accounts  for 

-  Ambiguities  in  MFP  output 

-  Signal  model  mismatch 

-  Colored  noise  (discrete  interferers) 

-  Adaptive  training  effects 

•  Challenge:  previously  no  way  to  achieve 
above  for  specific  adaptive  beamformers 
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•  Boundary  for  acceptable  MFP  performance  often  occurs  in 
threshold  region  =>  must  predict  threshold  region  accurately 
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